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Observational evidence exists that winds of massive stars are clumped. Many massive star 
systems are known as non-thermal particle production sites, as indicated by their synchrotron 
emission in the radio band. As a consequence they are also considered as candidate sites 
for non-thermal high-energy photon production up to gamma-ray energies. The present work 
considers the effects of wind dumpiness expected on the emitting rclativistic particle spectrum 
in colliding wind systems, built up from the pool of thermal wind particles through diffusive 
particle acceleration, and taking into account inverse Compton and synchrotron losses. In 
comparison to a homogeneous wind, a clumpy wind causes flux variations of the emitting 
particle spectrum when the clump enters the wind collision region. It is found that the spectral 
features associated with this variability moves temporally from low to high energy bands with 
the time shift between any two spectral bands being dependent on clump size, filling factor, 
and the energy-dependence of particle energy gains and losses. 



1 Introduction 

Evidence for particle acceleration to relativistic en- 
ergies mediated by the supersonic winds of massive, 
hot stars comes from the observation of non-thermal 
radio emission (e.g., Abbott et al. 1986). This has 
been interpreted by synchrotron emission on the 
basis of the measured spectra (much steeper than 
the canonical value ar ^ -1-0.6, F^, oc i'"'') and 
high brightness temperatures of ~ 10^~''K, far ex- 
ceeding lO'^K expected from free-free emission 
from a steady-state isothermal radially symmetric 
wind (Wright & Barlow 1975). Those particles has 
been suggested to be accelerated either in shocks 
caused by the instability of radiativcly driven winds 
(White 1985), in the shocked wind collision region 
of multiple systems (e.g., Eichler & Usov 1993) or 
in the termination shock (Volk & Forman 1982). 

Through a statistical study the presence of 
non-thermal radio emission has been linked 
to the binarity status of the stellar systems 
(Dougherty & Williams 2000), which is in support 
of a scenario where particles being predominantly 
accelerated at the forward and reverse shocks from 
the colliding supersonic winds from massive stars. 
Therefore, in this work I am considering the non- 
thermal particle spectra from being produced in the 
collision region of a typical long-period massive bi- 
nary system. 

By now ample evidence has accumulated that hint 
towards the dumpiness of WR- and 0-star winds: 
non-converging mass loss rates from various meth- 
ods (thermal radio. Ha, UV- lines. X-ray and IR di- 
agnostics, etc.), polarization and photometric vari- 



ability, stochastically variable substructures in lines 
(e.g., Eversberg et al. 1998), the observation of non- 
thermal radio emission even in some short period bi- 
naries (e.g., CygOB2#8A: Blomme 2005) where self- 
absorption in the optically thick winds should pre- 
vent the visibility of a non-thermal component, etc. 
Still open, however, remains a qualitative assessment 
about the filling factors, typical clump sizes, clump- 
ing startification, etc. 



2 Emitting electron spectrum 

This work considers, for the first time, the effect of 
dumpiness in the colliding winds of binary systems 
on the resulting non-thermal electron spectrum in 
the collision region at a given orbital phase, assum- 
ing that the particles to be accelerated stem from 
the pool of thermal wind particles. It requires a 
fully time-dependent treatment of particle injection, 
acceleration and losses. 

For this purpose I describe diffusive particle accel- 
eration and losses (radiatively, and via escape from 
the collision region with rate Tq^) that governs the 
non-thermal emitting electron spectrum, by the ki- 
netic equation: 

|(iV(i?,0) + ^{EN{E,t)) + = Q{E,t) 

(1) 

where E = aE — E\oss , the radiative energy loss rate 
E\oss = (^syn +^'ic)^'^ from synchrotron and inverse 
Compton losses in the Thomson regime, the acceler- 
ation rate a = ^"3^'^^^^"' (Reimer et al. 2006) with 
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Cr the shock compression ratio, and Ka the diffu- 
sion coefficient, Tq = ^ the escape time scale for 
a constant post-shock flow velocity V, and the size 
of the acceleration region, assumed cylinder-shaped 
for simplicity, with vq = K^jV . In this picture the 
dumpiness of the wind translates into a variable in- 
jection rate Q{E,t). while for a homogcnuous wind 
Q{E,t)oz6{E~Eo). 

In the following I consider the typical setting of a 
colliding wind region (see e.g., Eichler & Usov 1993) 
at a given orbital phase with a clumpy wind of vol- 
ume filling factor /vol defined as the total volume oc- 
cupied by the clumps with respect to the total wind 
volume, /vol = V^iump/^ot- The simplified picture of 
a constant particle density within the clump and no 
particles in the interclump medium used here is suf- 
ficient to expose the basic properties of variability in 
a non-thermal component from a clumpy wind. Fur- 
thermore, /vol and VoB shall be constant sufficiently 
close to the shock region. 

Q{E,t) = Qo6{E - Eo)H{t - to)H{U - t) de- 
scribes the particle injection into the shock region for 
a clump of size ^ciump = VbB(i« — io)- The solution is 
analytic, and shown in Fig. 1 for tu — > cxo. With the 
spectral index s of the emitting electron spectrum, 
oc E~^ , determined by the acceleration rate and es- 
cape, the exponential built-up of the electron spec- 
trum up to a maximum energy Ec = a/ {bgyn + he) 
(see Fig. 2, soHd Hue) is governed by both, energy 
gains and losses. Steady-state is reached on hours 
time scale in the typical settings of the colliding wind 
region in long-period binary systems. 

The interclump phase is described by Q{E,t) oc 
E~'^6{t — tv) as the source term in Eq. 1. For a > 
in this phase the solution of Eq. 1 implies a contin- 
ued acceleration with the maximum electron energy 
increasing further exponentially (see Fig. 2, dashed 
line): the particle acceleration appears as a multi- 
(here: two)-stage process. For a = further accel- 
eration ceases and the emitting electron spectrum 
declines on typically hour times scales in the astro- 
physical environment considered here. Variations of 
the acceleration rate a may be plausible considering 
that possible variations of the shock conditions (e.g. 
compression ratio, etc.) as the clumpy winds collide 
re not unexpected. 

Combining injection and interclump phases and 
taking into account past clump injections into the 
shock results in a temporally evolving electron spec- 
trum in the wind collision region. An example for 
^ciump = lO^^cm is shown in Fig. 3. Clearly visible 
are the spectral features from injection phases, and 
they propagate from low to high energies as time 
progresses. Thus the flux variations from clumpy 
particle injections into the acceleration process in 
the shock at low energies are anticipated to precede 
the flux variations at high energies. Ultimately, the 
underlying physical reason for this time shift is the 
unavoidable energy dependence of the radiative par- 
ticle energy losses, and their gains in the acceleration 



process. In addition to the clump size and volume 
filling factor, the expected time shift is therefore de- 
pendent on the electron energy losses and gains. 
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Figure 1: Injection phase: Resulting emitting elec- 
tron spectrum at times 1, 3, 5, 7, 10, 
15 and 20 ksec (correspondingly, lower 
to upper lines) after to = in re- 
sponse to injecting particles into the ac- 
celeration process starting at to- Pa- 
rametcrs used for all figures presented in 
this work: bolometric luminosity Lqb = 
IO^Lq and effective temperature Teff = 
45000K of the star closer to the stagna- 
tion point, terminal wind velocities Vwr = 
VoB = 2000km/s, mass loss rates of the 
OB, Mob = lO-^Mg/yr, and WR-star, 
-^WR = 10~^MQ/yr, constant magnetic 
field in the shock region B=l G, a stellar 
separation of D = 5 • lO^'^cm, postshock 
flow velocity V = O.GVob, diffusion coeffi- 
cient Ka = lO^^cm^/s, shock compression 
ratio Cr = 4 and Eq = IMeV. Steady-state 
is reached on hours time scale. 



3 Discussion and Outlook 

In this work consideration has been given to the ef- 
fects on the non-thermal electron spectra, as pro- 
duced in a diffusive shock acceleration process from 
the pool of thermal wind particles, when the mas- 
sive star winds are clumped as they collide. In this 
case, the electron injection into the shock varies with 
time, dependent on the volume filling factor and typ- 
ical clump size. Consequently, flux variations are ex- 
pected in non-thermal particle and photon spectra, 
in addition to the thermal emission component. I 
have shown that those flux variations are expected to 
appear time shifted when comparing various spectral 



2 



Clumping effects on non-thermal particle spectra 



bands, with the low energy flux variations preced- 
ing those at high energies. This clumping signature 
in non-thermal spectra probes wind clumping in the 
very vicinity of the collision region, and may possibly 
show up in multifrcqucncy observations at energy 
bands that arc dominated by non-thermal emission, 
provided sufficiently dense data sampling. 

A metallicity dependence of clumping in non- 
thermal spectra may be probed by comparing clump- 
ing signatures from non-thermal ion spectra, ap- 
pearing eventually in the gamma-ray domain only 
through TT^-decay photon emission, with those from 
non-thermal electron spectra. For this purpose we 
plan to apply the presented prescription to clumpy 
ion injections as well. 

The presented scenario describes a simplified 
clumping picture, which nonetheless proves to be 
sufficient to expose the basic properties of a vari- 
able particle injection into the acceleration process 
in colliding wind regions. This general prescription 
allows to also account for more complex clumping 
configurations (e.g. density structure within clump 
and/or interclump medium, irregular nested clump 
sizes on various scales, etc.), by modifying the injec- 
tion term Q{E,t) accordingly. 
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Figure 3: Temporal evolution of the total emitting 
electron spectrum in the wind collision re- 
gion with clumpy particle injection and 
clump size ^ciump = lO^^cm. A temporal 
step size of At = Iks is used. All spec- 
tra are separated artificially in normaliza- 
tion for optimal visibility, a =const at all 
phases is used here as an example. 
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Figure 2: Exponential increase of maximum energy 
as a function of time with I ~ V{t — to). 
During the clump injection phase the max- 
imum energy reaches a level -g — ^~f(t) ^ 

i - - (-«^^) (solid 

line). During the interclump phase 
the maximum particle energy increases 
further if a > as — ^ = 

(dashed line). £'max,i = 10^ MeV and the 
rate of energy gain a =const at all phases 
is used here for demonstration. 
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